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We derive a straightforward variational method to construct embedded soliton solutions of the 
third-order nonlinear Schodinger equation and analytically demonstrate that these solitons exist as 
a continuous family. We argue that a particular embedded soliton when perturbed may always relax 
to the adjacent one so as to make it fully stable. 
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The distortion of pulses propagating through waveguides arises due to dispersive spreading. However, if the waveg- 
uides is made using materials the refractive index of which depends nonlinearly on the intensity of the propagating 
pulses, it will be possible to compensate the effect of dispersion by a pulse-narrowing effect [lj. A typical example 
, in respect of this is provided by optical fibers. Here the weak nonlinearity of the index of refraction arising due to 
t-H ■ Kerr effect produces a self modulation which in turn causes steepening of the wave. The dynamical interplay between 
the dispersive and narrowing effects produces the so called optical solitons. For an isotropic medium with cubic 
| l 1 nonlinearity the complex optical wave field <f>(x, t) is governed by the nonlinear Schrodinger equation (NLS) 

GO 

ifa + fat + H 2 ^ = o, 4> = <f>{x,t). (1) 

(-h ' Here t stands for time and x, the coordinate along the fiber. The suffixes x and t of 4> denote partial derivatives with 

l— J. d 2 4> 

respect to these variables. In particular, fat = ~q~2- The main reason for the robustness of the solitons described by 

(1) is that the wave numbers of the solitons lie in a range that is forbidden for the linear dispersive waves Q- As a 
result, linear wave can not be in resonance with the soliton to receive energy from the latter. We then have a stable 
soliton. 

The propagation of picosecond pulses is well described by the NLS equation in (1) which accounts for only the 
^\ ' second-order dispersion and self modulation. But for femtosecond pulses another important physical effect namely, 
the third-order dispersion comes into play. In this case the appropriate evolution equation for the pulse propagation 
is given by 

ifa + fat + \<P\ 2 fa= i(3fat- (2) 

* ^ | Understandably, the parameter (3 is a measure of the perturbation caused by the third-order dispersion. Due to the 
effect of perturbation a soliton could have wave numbers residing inside the linear spectrum of the system. Such a 
soliton is called an embedded soliton £|. Clearly, the embedded solitons can be in resonance with the perturbation 
to generate continuous wave radiation for subsequent decay. Despite this, in certain parameter regimes these solitons 
can be stable [5|- In view of this the stability analysis of embedded solitons arising in the solutions of higher-order 
NLS equation has been subject of considerable current interest. The present work is an effort in this direction. 

Studies in embedded solitons [|| within the framework of a third-order NLS equation indicate that these soliton exits 
as a continuous family. The members of the family being characterized by their propagation velocities are connected 
through a gauge transformation 0. We demonstrate the existance of such a family by using a straightforward 
analytical model and argue that embedded solitons can be fully stable despite their tendency to shed radiation. We 
believe that such an approach has distinct advantages to deal with physical problems because many unknown effects 
are then readily expressed and evaluated. To achieve our goal we shall derive a variational approach involving trial 
functions in order to describe the main characteristics of the embedded solitons as determined by the third-order 
NLS equation in (2). Although the chosen trial functions have a specific form , the shape parameters are allowed 
to evolve as the solitons propagate. With this assumtion the evolution equation under consideration simplifies to a 
reduced Lagrangian problem \j\ . We shall try to study the stability of embedded solitons by taking variations of this 
Lagrangian with respect to the shape parameters. 
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II. VARIATIONAL DESCRIPTION OF THE THIRD-ORDER NLS EQUATION 

For the NLS equation in (1) there exists a well-defined spectral problem Q such that one can write a closed form 
analytical solution of it in terms of sech functions. As opposed to this, the third-order NLS equation can not be 
solved analytically. Keeping this in view we look for a variational treatment of (2) to find an accurate approximation 
solution. For our treatment we will essentially rely on a Ritz optimization procedure applied to the Lagrangian 
function for the NLS equation in (2). We have found that the action functional 



W = 



with the Lagrangian density 



(3) 



C 



(4) 



reproduces the third-order NLS equation via the Hamilton's variational principle. 

In the Ritz optimization procedure, the first variation of the action functional is made to vanish for a suitable 
chosen trial function. In analogy with the well-known solution of (1) we introduce a sech type trial function 



4>(x, t) = 7y(a;)sech 



t - y(x) 

a(x) 



exp i 



V(x)(t - y{x)) + (t - y{x)f + a(x) 



(5) 



Here ry, y, V , a, b, and a are real functions of x. The parameters 77, y and a are related to the three lowest-order 
moments of the <fi envelope and represent respectively its amplitude, central position and width respectively. The 
other parameters cr, V and b stand for the phase, velocity (centre of the soliton) and frequency chirp. Understandably, 
these parameters will all vary with the distance of propagation. Using (5) in (4) we get 



£■8 ^ t £s 5 



(6) 



where 



(1 \ o I dV 1 1 — y\ dy if db da \ / 1 — y \ dy 1 1 — y\ da 
C s ~ V < a— ■ ' 



, . - V— + - a b— 

dx \ a I dx 2 V dx dx 



dx V a 



— \ sech 
dx 



21 t-y 



(7a) 



4 2) = -~V 4 sech* 



t-y 



(7b) 



rtanh 2 ' 



V 



b(t - y) 



sech 



2 1 l -y 



(7c) 



and 



44) = _ N 2 f ]!_ l(uih (t_y) l{ r . b ±lA \ I ' - 



a a* 



a z a 



sech 2 ' 



V 



Kt - y) 



}) 



sech 2 ' 



(7d) 



Here the subscript s on C indicates that we have inserted the sech ansatz for <j)(x, t) into the Lagrangian density. In 
terms of (6) the variational principle implied by (3) leads to 



5 / (L) dx = 



(8) 



3 



with 



We have found that result for (L) is given by 



r+oo 

(L) = I C s dt. (9) 



.dy , tt 2 ( db u da\ x do 2 2 2 2 tt 2 6 2 / 2 2 , tt 2 



W=°^ - 2 ^ + 12 hfe" 6 ^ + *T3" 2+ 3^ + 2y 2 + ^T+^ + ^ ' + — ) I (10) 



III. PARAMETERS OF THE TRIAL FUNCTION 

From the vanishing conditions of the variationals ^j^-, ~52T' ~W~ an d ~W~ we obtain the following 

equations. 

Tx ^ - °' (lla) 

„ T dy tt 2 / d& ,da\ dcr 4 , 2 „ tt 2 6 2 mr ( 2 - tt 2 6 2 \ „ 

^(WV)=0, (lie) 



2t/ dy tt 2 / d& & da\ + do 2 2 2 2 ^ 2 tt 2 6 2 
dx 12 I da; da;/ da; 3^ 3a 2 6 



and 



- 2U - /3 — + 3U 2 + — 
dx V a 2 4 



+ w + — =o. (iif) 



From (11a) 

Use of (12a) in (11c) gives 



arj 2 — constant = E (say). (12a) 
V = constant. (12b) 



Again (12a) and (lid) give 

dy ir 2 a db Tr 2 b da da 2 , 2 , 9 7r 2 6 2 „ ( 9 7r 2 6 2 2 . 

- 2V I + -Tx-l2Tx + Tx-^ -r +2V + — + P V { 2V + — -^)=° ^ 

Similarly from (12a) and (lie) we have 



^ = rb, r = 2 + 6(3V. (12d) 

ace 
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Subtracting (116) from (12c) we write 

db 4 ( r A n<n 

T^^-aK-i- 11 )- (13) 

From (12a), (12d) and (13) we get a second-order differtial equation for the pulse width 



d 2 a _ 4r / r E 
dx 2 7T 2 \ a 3 a 2 



(14) 



The first integral of (14) gives 



where the potential field 



1 /rfa N 2 



2 V dx 



11(a) = (15) 



n(a) = ^-^ + K. (16) 
n z a z n z a 

The constant of integration K is determined by the initial conditions of (14). 

IV. POTENTIAL FUNCTION FORMULATION FOR THE DYNAMICS OF EMBEDDED SOLITONS 

First from the initial conditions of (14) written as a(x)\ x =o — a o an d (^|) | x= o = 0, we have 

K = -^ + ^. (17) 

Using (17) in (16) we write 

2r 2 ( ai \ 4rE /a 



n(a) = 4»H"l -^(--0- (18) 
Tr z a z \a z ) Tr z a V a / 

We now introduce a normalized pulse width z(x) = ^) to rewrite (15) as 



1 / dz s 



with 



2U +n w = (19) 



U{z) = (1 + 3/3V) I - - l) lp (1 + 3/3V) f - + 1 ) + v \ (20) 



z j \ \z 

and 

8 8En 



2 4 2 2' 



(21) 



The potential function H(z) vanishes for two values of z, namely, z\ = \ and z 2 = — ^+^+3)^/) implying that there 
is a minimum in between these points. The point at which H(z) is minimum can be obtained from dT ^^ = as 

z rn = . (22) 

In the limit when - = —2 (1 + 3/3V) we have z^ = 1 = z m and ^^|z=i = n(z m ) = 0. In this case the potential well 
degenerates into a single point such that a particle released at that point will stay there. In the context of this work, 
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FIG. 1: 11(2;) as a function of z for - 



2 (1 + 3/3F) 



this signifies that a wave pulse for which — = — 2 (1 + 3/31/) propagates with unchanged shape as a cosequence of exact 
balance between nonlinear and dispersive effects. Interestingly, we now have a continuous family of soliton solutions 
because - depends on the continuous variable [3V . This implies that in the neighbourhood of each embedded soliton 
there is another one having a slightly higher or lower energy. A particular embedded soliton when perturbed may 
always relax to the adjacent one. This makes an embedded soliton fully stable. 

In Fig. I we plot the potential function H(z) as a function of z for fiV = —0.2,-0.1,0,0.1,0.2. Each of these 
potentials has z\ — Z2 = z m = 1 and is, therefore, associated with an embedded soliton. The solid curve (3V = is 
the potential function corresponding to the usual second-order NLS equation while the dotted curves represent the 
potential functions of the embedded solitons as found in the solution of third-order NLS equation. For /3V > 0, the 
potential curves for the embedded solitons lie above the solid curve for z > 1. For (3V < the dotted curves fall 
below the solid curve. We have drawn only a few curves for discrete values of (3V . However, a dense set of curves can 
be drwan by varying (3V continuously. 
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